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UNCLASSIFIED  ABSTRACT 


EQUILIBRIUM  STRUCTURE  OF  THIN  TR-669(S6240- 10)- 1 

DIFFUSION  FLAME  ZONE,  by  P.  M.  Chung  December  1965 

and  V.  D.  Blankenship 


The  method  of  inner-and-outer  expansions  was  used  to  analyze  the  initial 
broadening  of  a  diffusion  flame  sheet  caused  by  the  nonvanishing  equilibrium 
constant.  Analytic  expressions  are  derived  for  the  thickness  of  the  flame 
zone,  and  the  position  and  the  value  of  the  maximum  flame  temperature.  The 
thickness  of  the  flame  zone  created  by  the  nonvanishing  equilibrium  constant 
is  found  to  be  less  than  that  created  by  the  decreasing  Damkohler  number 
studied  previously.  Varying  the  stoichiometric  coefficient  of  the  fuel  is  found 
to  affect  the  oxidant  side  of  the  flow  to  a  greater  degree  than  the  fuel  side 
itself,  and  vice  versa.  The  problem  of  nonvanishing  equilibrium  constant  is 
important  when  a  fuel  is  injected  into  a  hypersonic  stream  and  the  high  energy 
of  the  stream  begins  to  prevent  the  full  release  of  the  heat  of  combustion. 
(Unclassified  Report) 
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I.  INTRODUCTION 


In  the  combustion  of  initially  unmixed  reactants,  when  it  is  considered 
that  the  chemical  kinetics  is  represented  by  the  simple  relationship 

kf 

a  (oxidant)  4-  b  (fuel)  ^  d  (product),  (1) 


*b 


it  is  well  known  *  that  a  diffusion 
rf  =  Damkohler  number  = 


flame  sheet  is  established  as 

characteristic  residence  time 
characteristic  reaction  time 


oo.  (2) 


and  at  the  same  time 


equilibrium  constant 


0. 


For  this  limiting  case,  the  combustion  is  completely  controlled  by  the  diffusion 
and  heat  transfer  characteristics  of  the  flow  field  alone.  This  independence  of 
the  combustion  on  the  chemical  reaction  kinetics,  together  with  the  fact  that 
the  equilibrium  constant  is  zero,  enables  one  to  postulate  a  certain  set  of 

conditions  at  the  flame  sheet,  and  to  solve  the  governing  equations  rather 

.1,2 

simply 

When  either  the  Damkohler  number  or  the  equilibrium  constant  is  finite 
and  nonvanishing,  a  zone  of  finite  thickness  instead  of  the  thin  sheet  is  required 
for  combustion.  Solution  of  the  governing  equations  is  then  no  longer  simple, 
and  usually  only  a  full  numerical  treatment  can  satisfactorily  describe  the 
burning  process.  Such  numerical  analyses  and  certain  empirical  simpli- 

3 

fications  derived  from  the  numerical  results  were  given  by  Chung,  et  al  and 

4 

Libby  and  Economos  . 
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The  most  efficient  combustion  from  the  standpoint  of  heat  release  takes 
place  as  r.  -*>  oo  and  k^/k^  0  and,  hence,  as  the  thin  flame  sheet  is 

established.  It  is,  therefore,  very  helpful  if  one  can  quickly  predict,  without 
having  to  solve  numerically,  the  values  of  the  Damkohler  number  and 
equilibrium  constant  at  which  the  flame  sheet  first  begins  to  broaden  and  the 
flame  temperature  begins  to  drop.  Also,  one  would  gain  a  considerable  in¬ 
sight  into  the  development  of  diffusion  flames  by  studying  in  detail  the  manner 
in  which  the  flame  sheet  first  broadens,  and  by  studying  the  structure  of  the 
thin  flame  zone  thus  created.  The  above  information,  which  is  often  difficult 
to  obtain  through  a  pure  numerical  analysis,  can  be  deduced  by  a  suitable 
perturbation  analysis. 

Mathematically,  the  flame  sheet  represents  a  singular  point  *  since  the 

second  derivatives  of  temperature  and  reactant  concentrations  increase  without 

limit  as  one  approaches  the  sheet.  (See  Figure  1.)  Therefore,  the  flame 

sheet  broadening  cannot  be  analyzed  by  the  regular  perturbation  method.  In 

order  to  perturb  about  the  singular  point,  the  technique  of  inner-and-outer 

3 

expansions  was  successfully  employed  by  Chung,  et  al  when  the  broadening 

was  caused  by  the  decreasing  Damkohler  number  for  the  general  combustion 

process  given  in  Eq.  (1).  A  similar  method  was  shown*  to  be  applicable  when 

the  flame  sheet  broadening  was  caused  by  the  nonvanishing  equilibrium  constant, 

3 

with  Tf  -*  oo,  but  only  for  a  =  b  =  d  =  1.  Unlike  the  flame  zone  created  by  the 
decreasing  relaxing  the  restriction  of  a  =  b  =  d  =  lto  arbitrary  stoichio¬ 
metric  coefficients  does  not  result  in  a  simple  extension  of  the  analysis  for 
a=b=d  =  l  if  the  broadening  is  caused  by  a  nonvanishing  equilibrium  constant. 
It  is  because  the  flame  sheet  broadening  by  the  decreasing  directly  influences 
the  inner  region  only  whereas  that  due  to  the  nonvanishing  equilibrium,  constant 
directly  affects  both  inner  and  outer  regions  and,  hence,  influences  the 
matching  process. 


In  the  present  paper,  the  initial  broadening  of  the  flame  sheet  caused  by 
the  nonvanishing  equilibrium  constant,  for  equilibrium  flows,  wi)l  be 
analyzed  by  the  method  of  inner-and-outer  expansions  for  arbitrary  stoichio¬ 
metric  coefficients.  From  the  analysis,  closed  form  expressions  will  be 
derived  for  the  important  variables  such  as  the  value  and  the  location  of  the 
maximum  temperature  of  the  flame  zone. 


II.  FORMULATION 


The  basic  chemical  features  of  the  present  analysis  are  independent  of 

any  particular  flow  geometry.  We  shall  hence  consider  the  simple  stagnation 

3 

mixing  layer  configuration  previously  employed  for  convenience.  A  jet  of 
pure  gaseous  fuel  is  introduced  into  the  oncoming  uniform  stream  of  pure 
oxidant  at  the  stagnation  region  of  a  two  dimensional  or  axisymmetric  blunt 
body.  A  laminar  mixing  layer  is  then  formed  as  the  two  streams  meet 
(see  Figure  2). 

The  general  governing  equations  for  this  problem  have  been  derived  and 

appropriately  transformed  for  Prandtl  and  Schmidt  numbers  of  unity  in  the 

3  1 

previous  paper  ,  from  which  we  write  when  p-  =  0, 


f'" 

+  ff" 

=  0 

(H  + 

mj)" 

+  f  (H 
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+  f  (H 
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€  H 
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(3) 

(4) 

(5) 


/  e  \  / ,  a  Ml  b  M2  \ 

\  H  J  y1  “  "M^"  ml  "  TTJ"  m2 j  - 
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(6) 


where  c  is  the  suitably  normalized  constant  portion  of  the  equilibrium 
constant  and  it  is  defined  as 


c  = 


1 

a,  b 
a  b 


M' 


a+b- 


.  v  a+b  , 

‘(*)  (*) 


l-(a+b) 


K 


EO 


(7) 


with  the  equilibrium  constant  expressed  as 


kb  _  K 
17  "  keo 


Tn  exp 


("It)  • 


(8) 


The  stoichiometric  coefficient  d  is  set  to  one  in  Eq.  (1)  without  loss  of 
generality  because  Eq.  (1)  can  be  always  divided  through  by  d.then  a  and 
b  can  be  redefined.  In  Eq.  (8),  Kr.  is  a  constant.  The  exponent  n  is  set 
to  one  for  Eq.  (6)  .  AE,  R,  P,  and  T  are  the  heat  of  reaction  per  unit  mass, 
gas  constant  per  unit  mass,  pressure,  and  the  temperature  respectively. 

0  =  (AE  Cp)/(RAh°)  in  Eq.  (6)  with  the  modified  heat  of  reaction  Ah°  defined  as, 


Ah 


a  M 

■R7 


1  h° 
"  hl 


b  M- 

_ t 

to. 


(9) 


where  h  is  the  heat  of  formation  per  unit  mass,  c  is  the  constant  pressure 

P 

specific  heat  of  the  gas  mixture  which  is  assumed  to  be  constant.  In 

Eqs.  (3)  through  (5),  f(q)  is  the  conventional  Blasius  stream  function  which 

is  defined  such  that  f'  =  —  where  (  )'  denotes  the  total  differentiation 

e 

with  respect  to  the  similarity  variable  r\  defined  as 

*11/2 


*1  = 


(1  +°-) 

L(Pe  H 


(air  / 


P  dy 


(10) 


T n  =  0  was  employed  in  Ref.  3. 
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In  the  above  equation,  tr  is  either  zero  or  one  depending  on  whether  the 
flow  is  two  dimensional  or  axisymmetric .  u  is  the  x-component  of  velocity 
as  shown  in  Figure  2.  The  subscripts  e  and  -e,  denote  the  edges  of  the 
mixing  layer  on  oxidant  and  fuel  sides  respectively,  and  o  denotes  the 
stagnation  point,  p  and  p.  are  the  density  and  viscosity  of  the  gas  mixture 
respectively. 

The  subscripts  1,  2,  and  3  represent  the  oxidant,  fuel,  and  the 
combustion  product  respectively.  H  and  m  are  the  normalized  temperature 
and  mass  fraction  respectively  defined  as 


where  M  and  C  represent  the  molecular  weight  and  the  mass  fraction 
respectively. 

The  boundary  conditions  for  the  governing  equations,  Eqs.  (3)  through 
(6),  are  as  follows: 


For  q  -*-oo, 

f‘  =  1. 


H  =  H  =  — E-  T  , 


6  Ah°  e 


m3 

ml  =  TW[  =  mle* 


m2  =  0, 


(12) 
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for  rj  -  go, 


fl  = 

u-e 

u 

e 

H  = 

H 

-e 

H 

e 

ml 

=  0, 

m2 

m3 

b  M-, 

2,  -e  * 


(13)t 


and  for  r\  =  0, 
f  (0)  =  0. 


(14) 


Though  the  equilibrium  equation,  Eq.  (6),  is  algebraic,  it  automatically 
satisfies  all  the  boundary  conditions  because  the  fluid  is  always  in  chemical 
equilibrium  at  both  e  and  -e  due  to  the  fact  that  only  one  reactant  is  present 
at  the  mixing  layer  edges. 

V 

Now,  we  have  completed  the  formulation  of  the  boundary  value  problem 
which  is  comprised  of  Eqs.  (3)  through  (6)  and  boundary  conditions  (12) 
through  (14).  Solution  of  this  boundary  value  problem  constitutes  the  next 
section. 

III.  ANALYSIS 


The  boundary  value  problem  comprised  of  Eqs.  (3)  through  (6)  and  Eqs.  (12) 

through  (14)  will  be  analyzed  in  this  section.  The  solution  of  momentum 

5  £ 

equation,  Eq.  (3),  is  available  elsewhere  ’  for  some  values  of  \.  In  the 
present  analysis,  Eq.  (3)  was  first  integrated  for  various  values  of  \. 


2 

T  The  particular  flow  geometry  demands  that  H  /H  =  \  .  This  point  was 

3  -  e  e 

explained  in  detail  in  a  previous  paper. 
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When  c  =  0,  the  boundary  value  problem  can  be  solved  readily  by 

replacing  Eq.  (6)  by  a  set  of  postulations  derived  from  physical  reasoning 

at  the  flame  sheet.  Such  flame  sheet  solutions  which  a  priori  admit  certain 

singularities  at  the  sheet  are  well  known  ’  ’  ’  .  When  e  >  0,  on  the  other 

hand,  Eq.  (6)  is  highly  nonlinear  though  algebraic,  and  the  problem  is  in 

general  amenable  only  to  numerical  analyses.  The  present  interest,  however, 

as  was  explained  in  Section  I,  is  with  the  onset  of  the  flame  sheet  broadening 

and  the  structure  of  the  thin  flame  zone,  thus  with  the  small  values  of  €. 

We,  therefore,  shall  obtain  a  perturbation  solution  of  the  boundary  value 

problem  posed  with  c  as  the  perturbing  parameter.  As  it  has  been  pointed 

out  previously,  c  =  0  is  a  singular  point  and,  hence,  the  regular  perturbation 

technique  cannot  be  applied  to  the  present  problem.  In  the  following,  the 

method  of  inner-and-outer  expansions  will  be  employed  and  closed  form 

expressions  which  describe  the  salient  features  of  the  flame  sheet  broadening 

will  be  derived.  For  a  general  method  of  inner-and-outer  expansions,  the 

7 

readers  are  referred  to  Van  Dyke  . 


We  first  assign  an  inner  region  about  the  singular  point  (the  flame  sheet) 
and  represent  it  by  (i).  The  remaining  portion  of  the  flow  field  is  then 
designated  as  the  outer  region.  We  specify  the  outer  region  on  the  oxidant 
side  as  the  upper  outer  region  (u)  and  the  fuel  side  the  lower  outer  region  (L). 
The  superscript  (  )#  will  be  used  to  specify  various  values  at  the  flame  sheet 
which  exists  in  the  limit  of  e  =0.  (See  Figure  1.  ) 

For  the  two  outer  regions,  we  express  the  dependent  variables  by  the 
following  regular  series, 

ml,u  <€’  ^  =  Yu0  W  +  «j  Yul  W  +  ^  Yu2  <*■>  +  ° 
m2,u  (€’  Tl)  =  Fu0  (T1)  +  *J  Ful  {T1)  +  *2J  Fu2  (T1)  +  °  (*3J)’ 

Hu  (€ '  Tl)  =  Hu0  +  c  j  Hul  W  +  €  2j  Hu2  {r])  +  °  (€  3j)’ 


(15) 


and 


mitL(«.n)  =  yl0  (ti) +  €k  yl1  (T1)  +  «2kYL2(T,)  +  o(€3k). 

m2  L  (€ .  n)  =  flo  (t))  +  c k  fl1  (n)  +  € 2k  fL2  (n)  +  o  (c 3k), 


W  2k  3k 

Hl  <«,  t,)  =  Hlo  (q)  +«  Hl1  (r|)  +  ,  (r,)  +  O  («  5K). 


L2 


(16) 


For  the  inner  region,  we  first  define  a  new  independent  variable  with  a 
stretching  factor  as 


r  =  H  I  =  i_ 

b  q  q 

The  dependent  variables  are  then  expressed  by  the  series, 
Hij, i  (« ,  ;)  =  c1  y2  (o  +  «2i  $2  (4)  +  0(c3i), 
m2  i  (c.  4)  =  <i  Fj  (4)+€2i  F2(&)  +  0(€3i), 

H.  (c,  &)  =  H*  +ci  Hl  (&)  +  €2i  H2  (4)  +  O  (e31). 


(17) 


(18) 


Now,  it  is  crucial,  in  order  for  the  inner-and-outer  expansion  method 

to  work  with  the  proposed  series  (15)  through  (18),  that  a  particular  set  of 

constants  j,  k,  q,  and  i  be  found  such  that  the  ensuing  solutions  of  the  three 

regions  will  match  in  a  consistent  manner.  The  perturbation  parameter 

appears  explicitly  only  in  Eq.  (6)  among  the  governing  equations.  Since  c 

appears  in  Eq.  (6)  as  a  linear  parameter,  one's  natural  tendency  would  be 

to  expand  the  dependent  variables  for  the  outer  regions,  which  are  regular 

regions,  into  integral  powers  of  €.  Indeed  such  expansions  produced  correct 

asymptotic  solutions  for  the  outer  regions,  which  could  be  matched  to  the 

inner  solution  consistently  to  all  order,  when  the  flame  sheet  broadening  was 

3 

caused  by  the  decreasing  Damkohler  number  for  all  values  of  a  and  b.  Also, 
it  was  the  case  with  the  flame  sheet  broadening  caused  by  a  nonvanishing 


8 


equilibrium  constant  when  a  =  b  =  1.  A  similar  attempt,  however,  was  found 
to  fail  in  the  present  analysis  with  arbitrary  a  and  b.  The  physical  reason 
why  j  and  k  should  not  necessarily  be  an  integer  will  be  discussed  later.  A 
careful  study  of  the  series,  in  the  light  of  the  governing  equations,  and  the 
nature  of  the  matching  required,  showed  that  the  present  method  of  solution 
would  work  only  if  we  let 

j  =  1/b, 

(19) 

k  =  1/a. 

For  the  inner  region,  in  consideration  of  the  fact  that  this  region  becomes 
finite  as  the  second  term  of  Eq.  (6)  which  represents  the  backward  reaction 
becomes  non-negligible,  and  also  the  subsequent  matching  in  mind,  we  let 

q  =  i  =  1  /(a  +  b).  (20) 

The  above  values  of  j,  k,  i,  and  q  will  be  shown  to  be  correct  later 
through  matching.  The  series  developed  will  now  be  used  to  analyze  the 
various  regions  in  the  following  subsections. 

Outer  Regions 

The  series  (15)  and  (16),  with  the  values  given  in  Eqs.  (19),  are  first 
substituted  into  the  governing  equations,  Eqs.  (4),  (5),  and  (6),  and  a  set  of 
various  order  perturbation  equations  are  generated  for  the  outer  regions  by 
collecting  the  terms  with  equal  powers  of  c  .  The  perturbation  equations  of 
l-^s.  (4)  and  (5)  are  readily  integrated  in  terms  of  the  Blasius  function  f,  and 
the  boundary  conditions  specified  by  Eqs.  (12)  and  (13)  are  applied  to  the 
zeroth  order  solutions  for  the  upper  and  lower  regions  respectively.  The 
usual  homogeneous  boundary  conditions  are  satisfied  by  the  higher  order 
solutions.  By  applying  the  boundary  conditions  at  r|  -*  oo  to  the  upper 
solutions  and  those  at  r\  -*  -oo  to  the  lower  solutions,  one  of  the  two  constants 
of  integration  of  each  order  perturbation  solution  of  Eqs.  (4)  and  (5)  is 


eliminated.  The  first  two  order  solutions  of  Eqs.  (4)  and  (5)  and  the 
corresponding  perturbed  equations  of  Eq.  (6)  thus  derived  for  the  outer 
regions  are  as  follows: 


M 


1VJ>  *3 

Huo  +  Yuo  *  + 

(21) 

Hu0  *  Fu0  *  C  <£'  -  »>  +  He’ 

(22) 

<YuO,a<FuO)b  *  »• 

(23) 

Hul  +  Yul  *  A2  <£'  -  »• 

(24) 

Hul  +  Ful  *  C2  <£’  *  £>’ 

(25) 

'Yuo'*  <Ful>b  *  <HuO>a+b  (>  'TIT-  Y"°)eXP  1 

(26) 

HL0  +  YL0  =  B  (£'  -  M  +  H.e, 

(27) 

M , 

HL0  +  FL0  "  D  <£'  '  »  +  H-e  +  FT^  • 

(28) 

<YL0>a  <FL0>b  =  *• 

(29) 

HL1  +  YL1  »  B2  (f  -  M. 

(30) 

HL1  +  FL1  =  D2  <£'  - 

(31) 

<YLi>a(FL0>b  ■  <«L0>a+b  (‘  -  TI7  flo) 

('*&} 

(32) 

In  the  above  equations,  A,  C,  A^,  B,  D,  B2,  and  D2  are  the  remaining 
constants  of  integration  which  will  be  determined  subsequently  through 
matching. 


As  mentioned  earlier,  a  flame  sheet  is  established  when  «  =  0.  The 

well  known  flame  sheet  solution  is  obtained  for  the  present  problem  by 

applying  certain  boundary  conditions  at  r\*  to  Eqs.  (21),  (22),  (27),  and 

(28),  and  by  thus  determining  the  constants,  A,  C,  B,  and  D.  The  flame 
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sheet  solution  w£.s  given  in  the  previous  paper  from  which  we  write  for 
the  subsequent  use, 


*  Hto  * 


=  H*  = 


m,  m_  +  H  m ,  +  H  m_ 

le  2,  - e  - e  le  e  2,  -e 

m,  +  m_ 
le  2, -e 


Y'K  = 


Yio  ■ 


=  F*, 


*  Fio  =  °' 


f  (1V»)  =  X  + 


mZ,  -e  t1  -  X> 

mle  +  m2.-e 


Inner  Region 

A  substitution  of  the  series,  Eqs.  (18),  with  the  aid  of  Eq.  (20),  into 
the  governing  equations,  Eqs.  (4),  (5),  and  (6),  and  the  collection  of  the 
terms  with  like  powers  of  «  results  in  a  set  of  perturbation  equations  for  the 
inner  region.  The  perturbation  equations  generated  from  Eqs.  (4)  and  (5) 
are  integrated  readily.  The  first  order  inner  solutions  are  given  by, 


=  “  i  *  »■ 


A  A 

Hj  +F,  =  Vi  +  6, 


(Yjf^j)1*  =  (H*)a+bexp  (-  jgp), 


where  a,  (3,  y,  and  6  are  the  constants  of  integration. 

With  the  use  of  the  lowest  order  expressions  derived  for  the  outer 
and  inner  regions,  the  dependent  variables  will  now  be  matched  between 
the  regions. 


Matching 


In  view  of  the  governing  equations,  Eqs.  (4),  (5),  and  (6),  we  shall 

ct  b  i 

match  the  functions  (H  +  m^)  (H  +  rr^),  an<*  (mj)  (na^)  »  between  the  two 

outer  regions  and  the  inner  region,  to  the  lowest  order.  First,  by  matching 

the  functions  (H  +  m^)  and  (H  +  m^),  all  the  unknown  constants  to  the  lowest 

order,  A,  C,  B,  D,  a,  |3,  y,  and  6  will  be  determined.  With  the  constants 

&  b 

determined,  it  will  be  shown  that  the  functions  (m^)  (nr^)  also  match 

consistently  thus  completing  the  solution. 

The  function  (H  +  m^)  is  first  matched  between  the  upper  outer  region 
and  the  inner  region  in  the  following  manner.  We  obtain  for  the  upper  outer 
region  with  the  aid  of  Eqs.  (15), 


H  +  m, 
u  lu 


1 

■  <HuO  +  Yu0>  +  <E  <H 


ul 


+  Yul>  +  °(‘F) 


(39) 


Eq.  (39)  becomes  as  q  -*  q#,  with  the  aid  of  Eqs.  (33)  and  (34) 


H  +  m.  =  H#  +  £  (H*  +  Y'  )*+  V(H''  +  Y”)# 

u  lu  '  uO  uO  2  '  uO  u(r 


+  O  (43)  +  o(eF). 


(40) 


In  terms  of  the  inner  variable,  Eq.  (40)  becomes  with  the  aid  of  Eqs.  (17) 
and  (21), 


H  +  m.  =  A 
u  lu 


M  1  1  /  2  \ 

(f*  -  1)  +  He  +  C  Af"*  +  O  (cF)  +  O  U®75  ). 

1  (41) 


1  2 

Here  either  €F  or  ca+^  maybe  greater  than  the  other. 

The  function  (H  +  m^)  for  the  inner  region  becomes,  with  the  aid  of 
Eqs.  (18)  and  (36), 


H.  +  m.  =  H# 
i  i 


(42) 
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Eq.  (42)  becomes  in  the  outer  variable  as 


H.  +  m.  =  H* 
i  1 


+  Q  +  0 


A  comparison  of  Eqs.  (41)  and  (43)  shows  that  a  consistent  matching  up  to 

/  1  \  /  2  \ 

0\€  J or  Ou  J,  whichever  is  greater,  will  be  accomplished  for  function 
(H  +  mj)  if  all  the  explicit  terms  of  Eq.  (41)  can  be  matched  with  those  of 
Eq.  (43).  Eq.  (42)  is  equivalent  to  Eq.  (43)  except  that  the  former  is  in  the 
inner  variable.  We  therefore  match  the  explicit  term  of  Eq.  (42)  with  those 
of  Eq.  (41),  and  we  obtain 


A  = 


H  +  — TT— 

e  a  Mj 


a  =  A  f"#  , 


=  0. 


In  a  similar  manner,  the  function  (H  +  m^)  is  matched  between  the  upper 
outer, and  the  ipner  regions.  The  functions  (H  +  m^)  and  (H  +  m^)  are 
also  rhatched  to  the  lowest  order,  in  the  same  way,  between  the  lower 
outer  'regions  and  the  inner  region.  From  these  matchings  the  following 


results  are  obtained: 


H  -  H* 


C  = 


y  = 


(He  -  H#)  f"* 


6=0, 


B  =  A, 


D  =  C. 


Now,  all  the  constants  of  integration  for  Eqs.  (21)  through  (23),  (27) 

through  (29)#  and  (36)  through  (38)  have  been  determined  by  matching 

functions  (H  +  m^)  and  (H  +  n^).  Notice,  however,  that  the  matching  of 

these  two  functions  could  have  been  just  as  well  accomplished  with  many 

other  values  of  j,  k,  i,  and  q,  instead  of  those  employed  in  the  present 

problem.  We  shall  now  employ  the  crucial  equation,  Eq.  (6),  and  match 
ci  b 

the  function  (m^)  (m^)  *  an<*  thereby  show  that  only  the  particular  values 

of  j,  k,  f ,  and  q  given  by  Eqs.  (19)  and  (20)  permit  the  function  (m^)a  (m^)*5 
to  match  consistently  between  the  two  outer  regions  and  the  inner  region. 


We  derive  for  the  upper  outer  region  with  the  aid  of  Eqs.  (6)  and  (15), 


<mlu>a  <m2u>b  ■  *  <Hu</+b 


(*-^r  Y“o)eitp(-^)+°('1+F)  <46> 


The  above  equation  becomes  as  q  -*  q#,  (£  -*  0), 


<mlu>a<m2u)b  *  € 


('he) 


+  e 


a  4-K  ^  ^  1 

•<Hao >  U7  <y;0>*  «*p 


<HuO>*  . 

F75  9 '  exp 


(h;0) 


(-  4  )♦  «>  (*  -  t3~ySc)- 
('  Hf;)+<a+b)(HaJb'1)*<H;o>*  (>  -  ^7  ySo)- 


exp 


('4) 


+  o(tc) 


(47) 
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sjsp 


St 


which  becomes  with  the  aid  of  Eqs.  (21)  through  (23),  as  q  ,  q*  (J  -  0), 

<ml«,W  =  '  [(H>'>a+b  “P  ('  Ro)J  +  <  5  exp  (.  £)  • 

f  a  M. 

[“  "MJ"  (A“c)  (H*)a+b  +  C6  (H*)a+b-2  +  C(a+b)(H*)a+b“1 


+  Q(i2) 


For  the  inner  region,  we  derive  with  the  aid  of  Eqs.  (6)  and  (18), 


<mli>V2i>b  *  * 


<H*>a+be*P(-^)  +  «1+^  [.^L  ?,(«*, 


M 

"  b  Kq  F1  <H*>  +  (a+b)  (H*)a+b-1  Aj  +  (H*)a+b“2  h2  0 


(■  w>) 


Now,  before  we  can  analyze  the  behavior  of  the  above  equation  as  L  - 
we  must  f,»d  the  behavior  of  functions  ?,  «),  (;,  and  fl  f„  la 

values  of  5.  These  are  obtained  in  the  following. 

The  first  order  inner  expressions.  Eqs.  (36),  (37).  and  (38),  are 
first  manipulated  to  give 

*  •—  -  r  1 

(F.)  a  +  (a  -  Y)  l  (F .  )a  .  (H*)a+b  exn  (.  if  .  n 


x>  ►<> 


Solution  of  above  algebraic  equation  is  not  usually  obtainable  in  closed  form 
for  arbitrary  values  of  a  and  b.  However,  for  the  purpose  at  hand  which  is 
to  derive  analytic  expressions  for  Yj,  and  Aj  for  large  positive  values 
of  £,  an  asymptotic  solution  of  Eq.  (50)  is  obtained  as, 


Fj  (&)  = 


(H#)a+b  exp 


f  ip) 


*1* 


/^\  * 


(H*)a+b  exp 


(•  M 


(a  -  v) 


+  O 


,  A  A 

The  asymptotic  expressions  for  functions  Yj  and  are  similarly  derived 


j  te)  =  (a  -  Y) 


,  (H*)a+b  exp  (-  , 

4  +  - - ttt—vv-  — /J  +  O 


\a  ~  y 


(C)  =  v  t  - 


(H*)a+b  exp  (-  ^ 


1^ 


(a  -  v) 


Now,  having  derived  the  expressions  for  functions  Yj,  and  Aj  for 

large  positive  values  of  £,  we  return  to  Eq.  (49). 
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] 


With  the  aid  o£-Eqs.  (51),  (52),  and  (53),  Eq.  (49)  becomes  as  4  —  00  * 


(m1.)a(m2i)b  =  €  (H*)a+b  exp  (-  +  c 


(-A)- 


(a+b)  (H*)a+b_  1  y  +  (H«)a+b'2  0  v  -  (H*)a+b 


(a  -  y)  C 


,„lU\  /u-^a+b-l  ,  /TT,vxa+b-2  Q  ,  ,u.,.xa+b  a  Mi  /tT.1#%a+b  bM2 

(a+b)  (H+)  +  (Hv)  0  +  (H*)  —  +  (H-) 


[<H*)‘tb.  ”p  <r  ro)]E  .  +  0 


(a  -  y)  V 


Above  equation  becomes  in  terms  of  the  outer  variable, 


(m1.)a(m2i)b  =  c  exp  (H*)a+b  +  (a+b)  (H*)a+b“ 1  y 

+  (H" >  0  y  -  <H*)a+b  (a  -  V)  ”T7”“  6 

3 

+  o(,1+^)  ,  O  (,‘+^).  (55; 

Now,  the  outer  solution  as  q  -♦  q*  will  be  matched  to  the  inner  solution 
as  4  "*■«>•  This  can  be  done  by  matching  Eq.  (48)  with  Eq.  (55).  With 
the  help  of  Eqs.  (44)  and  (45),  which  define  and  relate  various  constants, 
it  is  readily  seen  that  all  the  explicit  terms,  of  Eq.  (48)  match  with  those 

cl  h 

of  Eq.  (55).  A  similar  matching  can  be  established  between  (m.  _  )  (m_  ) 

a  b  IE  <^E 

and  (nij.)  (m2.)  • 

Now,  we  have  accomplished  the  complete  first  order  matching  between 
the  two  outer  regions  and  the  inner  region. 


i 


With  the  present  values  of  j,  k,  i,  and  q,  the  higher  order  matchings 

can  also  be  carried  out  but  with  increasing  algebraic  complexities.  Also, 

it  is  noted  here  that  any  other  values  of  j,  k,  i,  and  q  do  not  permit  the 

&  b 

matching  of  the  function  (m^)  (m^)  * 

Maximum  Flame  Zone  Temperature  and  the  Location 

From  Eqs.  (18), 

1  / 

H.  =  H*  +  Hj  +  O  \( 


and  from  Eq.  (37), 

Hj  =  V  &  “  (57) 

In  order  to  obtain  the  first  order  inner  functions,  H,  (£),  Y,  (£),  and 
A 

Fj  (C)i  for  all  values  of  £,  the  algebraic  equations,  Eqs.  (36),  (37),  and 
(38),  must  be  solved.  This  entails  the  general  solution!  of  a  nonlinear 
algebraic  equation,  Eq.  (50).  For  most  cases  of  combustion,  however, 
the  primary  interest  is  with  the  value  and  the  location  of  the  maximum 
flame  temperature  as  affected  by  the  flame  sheet  broadening.  These 
informations  can  be  obtained  without  the  general  solution  of  Eq.  (50).  In 
the  following,  analytical  expressions  will  be  derived  for  the  value  and  the 
location  of  the  maximum  temperature  of  the  broadening  flame  zone. 

The  maximum  value  of  H.  exists  where 


T  Note  that  only  the  asymptotic  solution  of  Eq.  (50)  for  large  values  of 
obtained  analytically,  was  needed  for  the  matching. 


First,  we  shall  evaluate  (F.)  ,  where  (  )  denotes  the  value  at  the  position 

'  I'm  m 

of  the  maximum  flame  zone  temperature.  Note  that  r|  =  ti*  when  c  =  0.  We 

dFi 

differentiate  Eq.  (50)  with  respect  to  £  and  substitute  y  for  — in  accordance 
with  Eq.  (59).  The  resulting  equation  is  then  readily  solved  for  (Fj)  as 


<Fl>m 


(a  -  V)  V  (J) 

a  +  V  (f) 


m 


(60) 


Eq.  (60)  is  now  substituted  back  into  Eq.  (50)  and  the  resulting  equation  is 


solved  for  C 
temperature, 


m 


Then,  we  obtain  for  the  location  of  the  maximum  flame  zone 


Eqs.  (60)  and  (61)  or  (62)  give 


(1-F)^  then  follows  from  Eqs.  (56),  (57),  (61),  (62),  and  (63),  with  some 
manipulation,  as 


The  properties  at  the  thin  flame  sheet,  (  )#,  are  given  in  Eqs.  (33)  through 

(35).  Now  having  derived  the  analytic  expressions  which  describe  the  behavior 
of  the  maximum  flame  temperature,  the  results  will  be  discussed  in  the 
following  section. 


IV.  DISCUSSION 


First,  it  is  seen  from  Eqs.  (17)  and  (20)  that  the  order  of  the  thickness 

1 

of  the  flame  zone  created  by  the  nonvanishing  equilibrium  constant  is  c 

3 

In  the  previous  work  ,  the  order  of  the  flame  zone  thickness  created  by 

1 

decreasing  Damkohler  number  was  (€  ^)a+k+1  where  «D  is  the  reciprocal  of 


20 


Damkohler  number.  Thus,  the  thickness  of  the  flame  zone  created  by  a  finite 
equilibrium  constant  is  less  than  that  created  by  a  decreasing  Damkohler 
number  when  a  and  b  are  order  1. 


3 

It  was  found  that  all  the  perturbation  terms  in  the  outer  regions  except 

the  zeroth  order  terms  were  found  to  be  zero  in  the  analysis  of  flame  sheet 

3 

broadening  due  to  a  decreasing  Damkohler  number.  Hence,  there  the  outer 
regions  were  unaffected  by  the  broadening  flame  sheet  to  all  orders.  This  is 
the  reason  the  expansions  of  the  outer  functions  were  not  influenced  by  the 
stoichiometric  coefficients,  and  were  in  integral  powers  of  for  all  values 
of  a  and  b.  In  the  present  problem,  on  the  other  hand,  it  was  seen  in 
Section  III  that  the  expansions  of  the  outer  functions  depended  directly  on  a 
and  b.  Also,  it  was  seen,  in  that  section,  that  the  higher  order  outer  terms 
were  not  zero  but  would  have  matched  consistently  to  the  appropriate  inner 
terms.  Thus,  the  effect  of  broadening  flame  sheet  is  felt  in  the  outer  regions, 
as  well  as  in  the  inner  region,  when  the  broadening  is  caused  by  a  nonvanishing 
equilibrium  constant.  We  also  see  from  Eqs.  (15),  (16),  and  (19)  that 
increasing  the  stoichiometric  coefficient  of  the  fuel  supplied  from  the  lower 
outer  region  enhances  the  effect  of  the  present  flame  sheet  broadening  on  the 
upper  outer  region,  rather  than  the  lower  region.  Similarly,  the  stoichiometric 
coefficient  of  oxidant  supplied  from  the  upper  region  affects  the  lower  region 
more  than  the  upper  region  itself.  Eq.  (62)  shows  the  location  of  the  maximum 
flame  temperature  of  a  broadened  zone  with  respect  to  that  of  the  flame  sheet. 


We  see  that  the  sign  of  the  quantity 


aM^M^ 


determines  the  direction  to  which  the  position  of  the  maximum  temperature 
begins  to  move  at  the  onset  of  the  flame  sheet  broadening.  It  is  seen  that 
when  and  H  =  Hg  the  position  of  the  maximum  temperature  is  un¬ 

changed  and  is  at  r|*.  It  is  interesting  to  note  that  the  position  of  the  maximum 
flame  temperature  tends  to  move  toward  the  lighter  reactant  (reactant  with  a 
smaller  molecular  weight)  as  the  flame  sheet  is  broadened.  Also,  it  would 
move  toward  higher  temperature  reactant  source.  Moreover,  increasing  the 
stoichiometric  coefficients  a  and  b  will  accentuate  the  temperature  effect, 


as  compared  to  the  molecular  weight  effect,  on  the  movement  of  the  position 

of  the  maximum  temperature.  One  can  readily  manipulate  the  results  given 

3 

in  the  previous  paper  and  show,  when  the  flame  zone  is  created  by  the 
decreasing  Damkohler  number,  that  the  position  of  the  maximum  flame 
temperature  is  governed  by  the  parameter, 


b  M,  (H*  ■  H  )  -  a  M.  (H*  -  H  ) 

l  *  C  A  C 


Therefore,  when  the  flame  zone  is  being  created  by  the  decreasing  Damkohler 

number,  the  position  of  the  maximum  flame  temperature  depends  on  (aM. )  and 

3  1 

(bM?)  rather  than  and  M^.  Also,  in  that  case  ,  the  position  is  determined 
for  a  given  value  of  the  above  parameter  through  another  numerically  computed 
function.  Hence,  the  sign  of  the  above  parameter  does  not  necessarily  reflect 
the  direction  to  which  the  maximum  temperature  moves.  2 

Finally,  the  maximum  flame  temperature  is  given  by  Eq.  (64)  to  0(t  ). 

Eq.  (64)  first  shows  the  expected  result  that  a  nonvanishing  equilibrium 
constant  reduces  the  maximum  flame  temperature.  For  those  values  of  0/H* 
which  would  make  the  exponential  term  the  predominant  term  of  Eq.  (64), 
the  amount  of  the  flame  temperature  drop  increases  as  coefficients  a  and  b 
are  increased.  Also,  in  that  case,  the  actual  order  of  the  flame  temperature 
drop  as  well  as  the  flame  zone  thickness  (see  Eq.  62)  is 


In  closing,  it  is  pointed  out  here  that  it  is  understood  in  the  evaluation 

of  £  and  H.  that  H#  is  greater  than  H  and  H  .  Otherwise,  the 
m  im  e  -e 

maximum  temperature  is  not  created  by  the  combustion  and,  therefore, 

& _  and  H.  would  have  no  meaning. 

m  im  ® 


V.  CONCLUDING  REMARKS 


The  initial  broadening  of  a  diffusion  flame  sheet  caused  by  an 
equilibrium  constant  becoming  finite  has  been  analyzed  by  the  method 
of  inner-and-outer  expansions.  The  Damkohler  number  is  assumed  to 
be  infinitely  large  and,  therefore,  the  fluid  is  in  local  chemical  equilibrium 
everywhere.  Analytic  expressions  have  been  developed  for  the  thickness 
of  the  flame  zone,  and  the  position  and  the  value  of  the  maximum  flame 
temperature. 

The  various  effects  of  the  stoichiometric  coefficients  on  the  flam^j  zone 

broadening  have  been  discussed  in  detail  in  the  preceding  section.  As  a 

whole,  it  was  found  that  varying  the  stoichometric  coefficients  affect  the 

flame  zone  and  the  outer  regions  quite  differently  when  the  flame  zone  is 

3 

created  by  a  nonvanishing  equilibrium  constant  compared  to  the  case  where 
the  zone  was  created  by  a  decreasing  Damkohler  number.  Also,  the 
stoichiometric  coefficient  of  the  fuel  affects  the  oxidant  side  more  than 
the  fuel  side,  and  vice  versa. 

Now,  the  problems  of  diffusion  flame  sheet  broadening  due  to  either 
the  decreasing  Damkohler  number  or  the  nonvanishing  equilibrium  constant 
having  been  analyzed,  the  next  problem  worth  considering  may  be  the  problem 
of  simultaneous  perturbation  wherein  both  the  Damkohler  number  and  the 
equilibrium  constant  become  perturbed  from  their  respective  limiting 
values. 


INNER  (I)  , 

THIN  FLAME  ZONE 

UPPER  OUTER  (u)  ^  LOWER  OUTER  (L 


OXIDANT 


0 

f  (o)  =  0 


Figure  1.  Sketch  of  Flame  Sheet  Broadening 
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